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Abstract 

A description and a detailed comparison of the Mo and Tsai and the Dubna radia- 
tive correction schemes is presented. Numerical comparisons made in the kinematical 
region of the NMC high energy deep inelastic electroproduction experiment are dis- 
cussed. An overall agreement between the two approaches in the region of low x and 
high y, where the radiative corrections are largest, is better than 2%. 



1 Introduction 



It is well known that information on the nucleon internal structure is contained in the 
electromagnetic structure functions or in the one photon exchange cross section, cxi 7 , for the 
deep inelastic lepton scattering off nucleons. However, a determination of the one photon 
cross section from the data, which is a goal of electroproduction (i. e. electron- and muon 
scattering) experiments, demands excluding contributions from other electroweak processes. 
These processes account for a large fraction of the measured cross section, especially in the 
low x and high y region. They cannot be discarded from the measured cross section on the 
event-by-event basis; the measured differential cross section can instead be multiplied by a 
correction factor calculated theoretically. This is called a radiative correction procedure. 

The first radiative correction scheme was created in the sixties by L.W. Mo and Y.S. Tsai 
HI, H, |3| (MT scheme) in connection with the early SLAC electron scattering experiments. 
Another approach, originally formulated for the CERN muon scattering experiment planned 
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by the BCDMS Collaboration, was suggested by the Dubna group in the seventies [Q and 
upgraded later [||](D scheme). In the analysis of the data from the deep inelastic electropro- 
duction experiments, both schemes were extensively used (cf. e. g. ||), the MT scheme often 
in the exact (i. e. no 'peaking approximation') and upgraded version. However, the precision 
of the recent experiments is so high, see e.g.0, that also the radiative correction procedure 
has to be based on more precise theoretical calculations. Therefore understanding similari- 
ties and differences in the two radiatiative correction schemes is of ultimate importance for 
concluding about consistency of results coming from different experiments. 

The goal of this paper is to compare analytically as well as numerically the most upgraded 
versions of the two approaches. Our experience with the radiative correction procedure in the 
deep inelastic experiments carried out at CERN by the EMC, NMC and BCDMS Collabora- 
tions allows to point out the problems encountered during the application of this procedure 
in the data analysis. The previous comparisons between the two considered methods were 
limited to analysing the numerical results of the early versions of the schemes || or only 
dealt with a subset of the radiative processes ||. 

The paper is organized as follows. In Section 2 the deep inelastic kinematics and cross 
sections are defined. Sections 3 and 4, together with Appendices A and B contain a descrip- 
tion of the MT and D schemes respectively. A complete set of formulae is given for each 
scheme using the original notation. Only minor simplifications and changes are introduced 
for the sake of clarity. Useful relations between the respective notations are given in the Ap- 
pendix C. The schemes are then compared in Sections 5 (theoretical ideas) and 6 (numerical 
results) and finally a summary is given in Section 7. 



2 Basic Definitions and Kinematics of Deep Inelastic 
Scattering 

As mentioned above, a goal of electroproduction experiments is to extract the differential 
cross section in the one photon exchange approximation (fig. la) from the data. This cross 
section can be expressed in the following way by the structure functions Fi(x,Q 2 ) and 
F 2 (x, Q 2 ) of the target: 



d 2 o"i 7 (x,Q 2 ) 4:71a 2 



dQ 2 dx Q 4 



Mxy\ F 2 (x,Q 2 ) ( 2m 2 \ 2 2 



In this equation a is the fine structure constant, m is the electron (muon) mass, Q 2 = —q 2 
where q 2 is the square of the four-momentum transfer between the incoming and outgoing 
lepton, x = Q 2 1 (2Mv) the Bjorken scaling variable, M is taken as the proton mass, E and v 
are the lepton's incident energy and energy transfer in the proton rest frame and y = v/E. 

The one photon exchange process described by the eq. (|1]) , is a part of the lowest order 
(or Born) electroproduction cross section, a B . The other part of a B proceeds from the Z° 
boson exchange. The two contributions cannot be separated experimentally. However in 
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the present fixed target experiments the involved virtualities (i. e. the Q 2 values) are small 
comparing to the Z° mass squared and therefore eq. ([!]) is a good approximation of the Born 
cross section. The Z° contribution will be discussed later in more detail. 

The differential cross section (1) can also be expressed in terms of structure functions 
F 2 (x,Q 2 ) and R(x,Q 2 ): 



d 2 a l7 (x,Q 2 ) 4vra 2 F 2 



dQ 2 dx Q 4 x 

where R is defined as: 



Mxy ( 2m 2 W(l + 4M 2 x 2 /Q 2 ) 
~ V ~^E~ + [ ~~Q 2 ) 2(1 + /?) 



(2) 



, 2l a L _{l + AM 2 x 2 /Q 2 )F 2 



R(x, Q 2 ) = — = > < - 1; (3) 

<jl and ot denote the cross sections for the longitudinally and transversally polarised virtual 
photon respectively. 

As mentioned earlier the radiative events account for a large fraction of the measured 
cross section, cr meas . These effects may lead to a wrong interpretation of the measured event 
kinematics. For example an elastic scattering from the target can be mistaken for a deep 
inelastic event if it is accompanied by an energetic bremsstrahlung photon not measured in 
the experiment. The magnitude of the radiative effects in the measured cross section will be 
characterized by the so called radiative correction factor, r)(x,y), defined as follows: 

r]{x,y) = —^. (4) 

^meas 

This factor is used in the data analysis |7j] and will also be employed in comparing results of 
calculations between different radiative correction schemes. 

Finally we have to stress that in the most of the deep inelastic experiments only the 
inclusive measurements are performed. This means that only incident and scattered leptons 
are measured and kinematics of the reaction is defined by the leptonic observables. How- 
ever the kinematics of radiative events cannot be defined by leptonic variables only, e. g. 
a bremsstrahlung photon emission by the lepton results in a measured Q 2 different from 
the actual one. Therefore, when appropriate, we shall make a clear distinction between the 
hadron- and lepton-defined variables. 



3 Mo and Tsai Scheme 



In the original MT scheme [IJ 0, Bl the following processes were considered to contribute 
to the measured electron-proton inelastic cross sectionF]: real photon bremsstrahlung from 

1 The original MT scheme was formulated for the electron scattering; therefore a large effort was put to 
quantify the effects of the energy loss in the electron passage through the target. We shall neglect them in 
this paper since in the examples of practical applications we shall deal with muon scattering only. 
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an initial and final electron (fig. lb), vertex correction (fig. lc) and vacuum polarisation 
correction (fig. Id). In the latter one only the electron and muon loops were originally 
considered. The proton structure was accounted for through the two unspecified structure 
functions. The MT scheme formulation is thus model independent. The mathematical 
formulation of the MT is non-covariant. 

Corrections depicted by the diagrams in fig. 1 can be divided into two groups: emission 
of real photons (fig. lb) with energy larger than A and those with energy smaller than A 
(fig. lb) together with virtual corrections (fig. lc,d). Contributions from the soft photon 
emission and from the vertex correction separately are infrared divergent but the divergences 
cancel when the contributions are considered jointly JI0" |. The parameter A may have a 
meaning of the energy resolution or in another words a maximal energy of the emitted 
photon which is still not detectable in the experiment. It is often called the 'infrared cut-off' 
parameter. This interpretation implies that A may not be too large. Too small A must also 
be avoided since it may cause numerical instabilities in the computation of the soft photon 
contribution. The numerical results of the calculations should not depend on A. 

In the MT approach the measured cross section can be expressed as follows (see Appendix 
A for the definitions of variables used in MT publications): 



d 2 a n 



dudQ 



+ 



dudQ dudQ 



(5) 



where 



(6) 



is a residuum of the cancellation of the infrared divergent terms and takes into account 
all soft photon emissions in the lowest order of a. It is a well known fact that the infrared 
divergencies cancel in each order in a JEDI and therefore it is possible to sum up contributions 
from all soft photon emissions ||11|| . The exponential factor, e~ SR ( A \ in eq.([5[) is a result of 
this summation. The function F(Q 2 ) contains all A— independent terms: 



F(Q 2 ) = 1 + ^ac + Co 



6 V 



8„ 



(7) 
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5 e '^ 

"vac 



^vtx 



-5 4m^ 
9 3Q 2 



1 



4 



$ cos - 



1 I 4 <M 

Q 2 



Q 2 



E s — E f 
E x 



'1 + AmlJQ 2 + 1 
j l + Am 2 e JQ 2 - 



(9) 
(10) 



2 given here is the full formula. The approximation 2a/n(— 5/9+l/3\nQ 2 /ml ^) holds for Q 2 ^> m 
(see e.g. eqs 5 and Al from ref. 0). For extremely low Q 2 formula (8) converges to Q 2 /15m 2 
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and $ the Spence function?]. The second term in eq.(|5|), cr tails , accounts for contribution from 
processes where the real photons of energy larger than A are emitted: 

details = dV(a; > A) = fM- dM dV ^ (n) 
dvdQ dvdtt JMj=M 3 dvdVl ' 



where M denotes a target mass, M™ ax = y M 2 — Q 2 + 2M(u — A) and a^ r is given by the 
formula (A. 24) in || and quoted in Appendix A; Mj is an effective mass of the hadronic final 
state. The integration over Mj in eq.(|TTD means that all the final hadronic states contribute 
to the cross section measured in a kinematical point (Q 2 , v\ fig. 2, so called 'radiative tails': 
elastic (Mj = M), resonance production (Mj = M res ) ... and deep inelastic tails, i.e. tails 
from the continuum. 

Equation (H) does not take into account double- (multi-) photon exchange reactions 
(fig. le) nor any radiative correction from the hadron current (fig. lf,g). These effects were 
estimated only for the elastic e — p interaction . Application of these calculations to the 
inelastic interaction is incorrect. Thus these results were discarded in practical applications 
of the MT scheme. Evaluation of these effects in the inelastic e—p interactions was not done 
in this scheme. To do this a model of the proton internal structure is necessary. The best 
framework presently is the quark parton model (QPM). 



In the upgraded version of the MT scheme, used in the analysis of the NMC results |L2 
the t + t~ and qq loops in the vacuum polarisation |TU| , |T5| ] and the virtual photon — Z° boson 
interference [Hi were also included. 



4 Dubna Scheme 



In the Dubna scheme the calculations of the deep inelastic processes are based on a mixed 
approach which uses both a model independent and the quark parton model treatment 
of the radiative corrections. The radiative corrections to the leptonic current (fig. lb,c) are 
calculated within the same model independent approach as in the Mo and Tsai scheme except 
the way of treating the soft bremsstrahlung photons. All other corrections are calculated 
using the QPM approach. 

In the inclusive electroproduction experiments the radiative and non-radiative events 
cannot be distinguished and therefore the cr meas should not depend on any parameter A, 
a property, which implicitly holds in Mo and Tsai scheme. The D scheme is explicitly A 
independent as a result of an integration over the whole bremsstrahlung photon phase space. 
The relevant procedure is described in detail in [14[ . 



Applying the QPM and the fact that quarks are point-like objects permit in principle to 
calculate all QED processes to all orders of a. This means that in addition to the diagrams 



3 In the original MT formulation |3| the l q 2 ' symbol was used in eqs @ - (0)- This variable, however, was 
then redefined to become the measured four momentum transfer, q 2 = (s — pf, and thus coincides with our 
definition of — Q 2 . The latter was thus used here for clarity. 
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in fig. l.a-d also the hadron current corrections (fig. l.e-g) can be evaluated. The diagrams 
of fig. l.e-g thus become as in fig. 3a-c. All these processes were taken into account as well 
as the lowest order weak corrections, not shown here, which can also be calculated in the 
framework of the quark-parton model. The outline of the radiative corrections calculations 



in the D scheme is given below; the details are given in refs |14J jig] . 
The measured cross section is now expressed as follows^]: 

dVrneas dV 



dQ 2 dx dQ 2 dx 

dV in . tail d 2 a IR 



^ dQ 2 dx dQ 2 dx 

+ ^ E E E c > K ( B ,P)[ V (B,p)Rl(B)+pA(B,p)Rt(B) 

^ B=-y,I,Zb=i,q Q,Q 
a2 

u u el. tails /, cy\ 

+ ~d&dx~- [ ] 

The first two rows of this formula represent the results of the model independent calculations 
of the radiative corrections to the leptonic current. The third row represents the quark parton 
model description of the lepton-hadron interactions (fig. 3a) and of the radiative corrections 
to the hadronic current (fig. 3b, c) as well as certain interference terms. Finally the last term 
in eq. ([12]), d 2 <T e i.taiis/dQ 2 dx, describes the elastic and resonance radiative 'tails'. 

The d 2 a B ' /dQ 2 dx cross section in eq. ( |T2"D denotes the full Born cross section for the deep 



inelastic scattering, i. e. the cross section containing both the one photon- and one Z° boson 
exchange contributions: 

d 2 <r B 27ia 2 y * 2 1 B 2 

d^~^' Q ¥' M) ' (13) 
where the functions S R {y,Q 2 ) are: 



S?(y,Q 2 ) 
Si(y,Q 2 ) 



Q 2 - 2m 2 , 

2[(1 -y)S 2 -M 2 Q 2 ] 
2Q 2 S(2-y). 



(14) 



with Q 2 = Sxy and S = 2ME. The functions Ai are given in Appendix B. 



The Sr in eq. (|12|) is responsible for those parts of the soft and hard collinear photon emis- 
sions which could be resummed to all orders using the covariant exponentiation procedure, 



11, 16 . It reads: 




1 In 



y 2 (i 



X) 



(1 - yx)(l - y(l - x))' 



(15) 



Additional factor a/n in the third row of eq. (|12[) is hidden in the definitions of the R functions. 
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The 5 VR correction factor in eq. flT2| ) is a remnant of the exponentiation and of the subtraction 
procedure used to disentangle the infrared divergent terms from the d 2 cr in . tai i/d(5 2 dx cross 
section |L4], see below. It thus contains the vertex correction, fig. l.c and is given by 



:VR 



'vtx 



a 
2^ 



In 2 



(1 - yx) 



(l-y(l-x)) 



y) 



(1 - yx)(l - y(l - x)) 



$(1). (16) 



The d 2 <7 in .taii/d(3 2 dx in eq. ( |i2|) describes the inelastic radiative tail for the lepton current 
correction only: 



d 2 CTin.tail 

dQ 2 dx 



2a 3 y 
Sx 



dQ 2 h dM 2 h J2Mxh,Q 



jSi{y,Q 2 ,yh,Ql) 



(17) 



In this formula is the invariant mass of the final hadronic system and the variables bearing 
a subscript ! /t' refer to virtual photon-target vertex in contrast to the variables measured in 
the inclusive electroproduction experiment. For their definitions see Appendix B where also 
the explicit expressions for the radiator functions Si are given. 

It is a well known fact that the d 2 <Ti n ta ii/dQ 2 da; cross section is infrared divergent. To 
regularize it a simple trick ('fixation procedure') of adding and subtracting an extra term 



-a 



d 2„B 



a 



dQ 2 dx dQ 2 dx 



dQ 2 h dM 2 h F IR (y,Q 2 ,y h ,Ql), 



to d cri n .taii/d<3 dx was employed [TT], [T3J]. In the added term an integration over a full 
photon phase space was carried out, resulting in the above given expressions for 5r and S VR . 
The subtracted term appears explicitly in the second row of eq.(|l^), so that the difference 
d 2 cTi n .taii/d(5 2 dx — d 2 o~ IR /dQ 2 dx in eq. (p~2|) is finite over the full kinematic domain of Q\ and 
Ml- This method is a key point of the D scheme, making it explicitely A independent. The 
function T IR is given in Appendix B. 



The third row in the eq. fll2| ) represents the quark parton model calculations and contains 
clearly visible vector (7 and Z°) and axial (only Z°) contributions. Index B runs over the 
photon exchange (7), the Z° boson exchange (Z) and their interference (/). Index b stands 
for the scattering with the single photon emission by the quark (q), fig. 3c, and its interference 
with the photon emission by the lepton (i). The double photon exchange (fig. 3a) as well as 
the vertex corrections on the quark line (fig. 3b) are also hidden there. Both scattering off 
quarks (Q) and antiquarks (Q) were considered 0. Coefficients q, are equal to Qq or Q^Qq 
for the process of bremsstrahlung photon emission from the quark or for the interference 
term between the photon emission from the lepton and from the quark respectively; and 
Qq are the charges of the lepton and of the quark given in the electron charge units. The 
sign factor p is defined as: p = p^Pq where p^Q = ± 1 for particle (antiparticle). For a 
detailed form of the coupling strength factors K(B,p), modified vector, V(B,p) and axial, 
A(B,p) couplings as well as of the functions R^' A (B) we refer the reader to ref. ]TB| . 

5 In the original formulation [|l5| the index b assumes also values and 1 which correspond to the contri- 
butions from processes in fig. la,b,c. Here they are included in the model independent parts of eq. ( |l2| ) 
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The problem of the soft photons' emission from the quarks can be uniquely solved since 
the quarks are not observed and thus the quark states are all summed up in the cross 
section. The infrared divergencies for initial- and final state quarks as well as the quark mass 
singularities for the final state quarks cancel in each order in a M, plj. Mass singularities 



associated with initial quark are included in the definition of the proton structure function. 

The last row in eq. d 2 cr e i. t aiis/d(5 2 dx, is in the D scheme treated essentially in the 

same way as in the MT except that it is formulated in a covariant way. Finally the vacuum 
polarisation, fig. Id, was taken into account via the 'running' a(Q 2 ) which in the Q 2 ^> m 2 
approximation (rrif stands for the lepton and quark effective masses), |l3j is defined as 
follows: 



«(Q 2 



a 



I E/ CfQjS 



f ■ 

vac 



where c/ and Qf are the colour factor and the electric charge of fermions / (c/ 
for leptons); '/' runs over all leptons and quarks. 



(19) 



We shall close this section with the following remarks: first, although not shown ex- 
plicitely in eq. (0), the weak loop correction contribution is also present in the D formu- 
lation. It is calculated within the QPM framework; details are given in Second, the 



0(a 2 ) corrections (a 4 contributions to the cross section) not shown explicitely, were also 
implemented. For the elastic radiative tail they were calculated completely in the first paper 
of ref. |J, while for the inelastic continuum they were implemented in an approximate way, 
described in the second paper of ref. 



5 Comparison of the MT and D Schemes 

Results of deep inelastic experiments were analysed using either MT or D radiative correction 
schemes. Therefore it is of ultimate importance to understand the differences and similarities 
between the two approaches. In this section we shall make a brief summary of theoretical 
ideas in the two schemes; numerical comparison will be presented in the next section. To 
faciliate the comparison, the relations between the variables used in the MT and D formulae 
are given in Appendix C. 

The D scheme is formulated in a covariant way and its model dependent part is based on 
the quark parton model. The covariant formulation means that all formulae are expressed in 
terms of the Lorentz invariants and are thus independent of the choice of the reference frame. 
This also means that the D scheme is explicitly independent of the infrared cut-off param- 
eter A, cf. eq. (p!2|). The A— independence was obtained through a special mathematical 
procedure. 

The MT scheme is not covariantly formulated. It should also be A— independent which 
formally means that the derivative over A of the right hand side of eq.(|^) should be equal 
to zero. The calculations show that this does not hold and that a certain dependence of the 
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results on A should be observed. For not too small A this dependence is very weak. For 
very small A (apart of numerical instabilities), eq.([5]) depends much stronger on the value 
of A which in this case should be calculated separately for every kinematical point (i.e. for 
every (x, y) value) in order to minimise the dependence. 

The elastic and resonance tails (i. e. the contribution of the reactions with the elastic and 
resonance final states) are calculated in the same way in both schemes. The inelastic tails 
(the contribution of the reactions with the inelastic final states) originating from the leptonic 
bremsstrahlung are treated differently but the differences are purely mathematical. The 
inelastic tail in the MT is calculated with the same formula as the elastic one (eq. (|TT|)) where 
the inelasticity of the process is taken into account through integration over kinematically 
allowed final state masses. The corresponding structure functions, Wl 2 {see Appendix A), 
are not specified; in practical applications the electromagnetic structure functions are used. 
The leptonic inelastic tail in D is also calculated with arbitrary structure functions while 
tails originating from the hadronic bremsstrahlung were calculated within the quark parton 
model framework. 

The usage of the parton model in the D scheme allows to calculate the hadron current 
corrections in the inelastic case as well as a double photon exchange process; they are cal- 
culated up to a 3 . The hadron current corrections cannot be calculated in the MT approach 
and therefore such corrections were not taken into account. The qq loop is also naturally 
present in the vacuum polarisation process in the D scheme while in the MT it was only 
added later. The same is true for the weak interactions contributions (Z° exchange and the 
7 — Z° interference). The weak contributions are however very small as compared to the 
present experimental resolution. 

Finally the a 4 lepton current corrections were partially taken into account in the D 
scheme but not in the MT. 



6 Numerical Calculations in the MT and D schemes 



A tremendous increase of the accuracy of deep inelastic electroproduction experiments de- 
mands a similar increase of the accuracy of the radiative corrections calculations. Therefore 
the early versions of these calculations have constantly been improved. At the same time 
comparisons between the two considered schemes were made, H, EL H. In the analysis of the 



recent very precise mesurements by the NMC ]7|, [12| both schemes were used in their most 
upgraded versions and applied for nucleon and nuclear targets. The MT code, named FER- 
RAD35, apart of the processes contained in the eqs (|5|)-(|rTD included also the tau lepton- 
and quark loops and the photon— Z°-boson interference. The tails were treated in an exact 
way in contrast to the peaking approximation 0, ||] applicable for the electron deep inelastic 
scattering, A detailed input information (structure functions, form factors, nuclear structure 
models, etc) was introduced |18|]. The D code, named TERAD86, was in principle used in 
the same version as sketched by the formulae ([12]) - (|19|) and employed in the dedicated 
BCDMS experiment H, except that the FERRAD35 input information was supplied. In 
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this section we shall compare the FERRAD35 and TERAD86 results on the radiative cor- 
rection factors, 7], as well as on the elastic radiative tails. Comparisons will be done in the 
kinematical region of the NMC positive muon-proton deep inelastic scattering experiment, 
i.e. for 0.003< x <0.9 and 0.1< y <0.9 at the incident muon energy 280 GeV. The input 
information in the calculations was as follows: Gari and Krumpelman proton form factor 
parametrisation, |19j structure function F 2 as measured and parametrised by the NMC, J/J, 
and finally for the R(x, Q 2 ) the parametrisation of SLAC was taken for all x and Q 2 >0.35 
GeV 2 . For smaller Q 2 the value of R was assumed to be constant and equal to the value at 
Q 2 = 0.35 GeV 2 . The proton resonances were neglected. For further details see refs [12j and 
0- 



The main numerical problem in the radiative correction programs is an integration of 
the radiative tails. For example, the integrand of eq. ([TT]) changes by 27 orders of magni- 
tude within the integration interval. It is obvious that this function should be integrated 
either by a high accuracy routines (e. g. the CERNLIB GAUSS routine demands an accuracy 
parameter e < 10 -12 ) or by dividing the integration interval into many small sections. It 
seems however that the best method would be to change the integration variable from cos 9k 
to log (— q 2 ) but also in this case high accuracy of the integrating routines is needed. This 
method was used in TERAD86 while in FERRAD35 dividing the cos 9k integration interval 
into many sections was normally employed (logarithmic integration in FERRAD35 was also 
tried and gave the same results). Due to the subtraction procedure in D the integrands in 
eqs (|T7|) and ([Tj|) are fairly smooth and do not demand any extreme precision of integration. 

The dependence of the FERRAD35 results on the parameter A was carefully studied. 
The results are presented in fig. 4. In the region of low x and large y where the radiative 
correction factor 77 is largest, the results are only weakly dependent on A for A >200 MeV. 
Results given below were obtained with A = 280 MeV. 

Comparison of the FERRAD35 and TERAD86 results is presented in figs 5-6. Radiative 
corrections are very large, exceeding 50 % at low x and high y, cf. fig. 5. In that region 
the agreement between the results of the two schemes is better than 2%, cf. fig. 6 (closed 
symbols). However the rf and qq loop contributions in the vacuum polarisation process 
(fig.l.d), absent in the original version of the MT scheme and included in FERRAD35, give 
up to 2% contribution to the radiative correction factor in most of the kinematic region (fig. 6, 
open symbols denote results of calculations without those contributions). Fluctuations visible 
in the high x part of the curves in fig.6 come from numerical instabilities of FERRAD35 in 
that region. 



7 Summary 

The two existing schemes of radiative correction procedure, the Mo and Tsai and the Dubna 
ones are differently formulated and are (partially) based on a different physics approaches. 
Both were extensively used in analysing the high energy experimental data. In this paper we 
presented the two schemes in detail and compared them analytically and numerically from 
the point of view of their effect on the results of the deep inelastic positive muon scattering 
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from a proton target at 280 GeV. To this aim we used the latest version of the D and the 
upgraded version of the MT programs. The latter included the r + r~ and qq loops in the 
vacuum polarisation and the virtual photon— Z° boson interference terms, all absent in the 
original formulation. In contrast to the 7 — Z° interference the quark loop contribution 
turned out to be quite substantial, changing the total radiative correction by about 2% in 
the measured region. 

The MT scheme contains the 'infrared cut-off' parameter, A. The results should not de- 
pend on its value (provided it is not too large and not too low) and indeed it is approximately 
so when A is equal to about 0.1% of the beam energy value. The covariant formulation of 
the D scheme excludes the existence of such parameter. 

The overall radiative correction reaches 50% at low x and high y. Calculated from the 
two schemes the corrections agree to better than 2% in this region. Differences are thus 
insignificant over the most of the phase space covered in the fixed target DIS experiments. 
They are of the order of other systematic errors in the data analysis [[/[] . 

Neither of the two radiative correction schemes contain a contribution from the multipho- 
ton exchange process to the elastic radiative tail which may be important for heavy nuclear 
targets. Results of quantitative estimates of those processes, relevant for the heavy target 
data currently analysed by the NMC, are discussed in a separate paper . 
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8 Appendix A 

Below we summarize variables used in the MT formulation. The metric used is such that 
ps = E P E S — ps and four vector components are in the laboratory system. Notation is 
explained in fig. 7 and the coordinate system is that of fig. 8. 



s = (E s ,s) four momentum of the incident lepton 

p = (E p ,p) four momentum of the scattered lepton 

9(fl) lepton scattering angle (solid angle), 

cos 9 = ~sp/ I s 1 1 p I 

9k , (pk bremsstrahlung photon emission angles 

9 v ,<fr v virtual photon emission angles 

t = (M, 0) four momentum of the target proton 
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k = (uj, k) four momentum of the bremsstrahlung photon 

Pf = s + t — p — k four momentum of the final hadronic system 

q 2 = (s — p — k) 2 = (pf — t) 2 four momentum transfer 

— Q 2 = (s — p) 2 measured four momentum transfer 



In the one photon exchange approximation and assuming one photon emission, the ra- 
diative tail from the j'th mass level can be written as (formula (A. 24) in @): 



d 2 aj tr 



where 



a 3 fE s 



2Mu;d(cos 6 k ) 
2n \E S J J -i q 4 (u — \ u \ cos9k) 



Wi(q 2 ) 



-am 



.r" 



2E s (E p + uu) + 



/ 2 

am 



2E P {E S -uj) + 



-2 + 2v(x~ 1 -y- 1 



m 



[s-p-to 2 ) + (s- p) [2E S E P - (s • p) + uj{E s - E l 



+x 



2 (E S E P + E s uj + E 2 )+^--(s-p)-m< 
-y~ l 2 (E S E P - E p u + E 2 ) + |- - (s ■ p) - m 2 j 
+W{(q 2 ) 

+4z/ (x" 1 - y -1 ) (s - p)(s-p - 2m 2 ) + (x~ l - y~^ (2s ■ p + 2m 2 - q' 



+ % J m 2 (2m 2 + q 2 ) +4 

X 6 yi 



UJ 

u 

u I 
u 2 

q 2 
a 

a 1 
b 
v 



- (u 2 - Mf) I (u - \u\cos9 k ) 



s + t — p 
E s + M- 



Pj + k 



■ E p 

1/2 



cos 6 P 
cos 6L 



x 



(2 2 

2m 2 + M 2 — 2(s -p) + 2M (E s - E p ) 

2m 2 — 2(s ■ p) — 2uj (E s — E p ) + 2oj \ u | cos Ok 

uj {E p — | p | cos # p cos 6^) 

uj (E s — | s | cos9 s cos 9 k) 

— uj I p I sin 6* p sin 9k 

(a' -a)" 1 

(| s | cos#— I p I)/ I w I 
(I s I — I p I COS 9)/ | m I 

,2\V2 



(a 2 -b 2 ) 



12 



y = (a'^-b 2 ) 1 



Wi, W 2 denote structure functions; in particular W( 2 (q 2 ) are the structure functions at four 
momentum transfer q 2 and invariant mass of the hadronic final state Mj. Wi t2 are connected 
with the functions F± }2 of eq. ([1]) in the following way: F 2 = vW 2 and F\ = MW\. Observe 
that the meaning of the u, x and y wariables used in the MT formulation is different from 
their generally accepted meaning as the DIS variables. 



9 Appendix B 



Below the exact expressions for certain functions in the D scheme will be given. The 'gen- 
eralized' structure functions Ai(xh, QfJ, in eq. ( |T3"D are: 



A 1 (x,Q 2 ) = 2F 1 NC (x,Q 2 ), 
A 2 (x,Q 2 ) = ±-F™{x,Q 2 ) 

yS 

A 3 (x,Q 2 ) = ^s NC &Q 2 ) 



with 



Fff&Q 2 ) = F 1 ^x,Q 2 ) + 2\Q e \(v e + \a e ) X (Q 2 )G 1)2 (x,Q 2 ) 

+ (v 2 e +a 2 e + 2\v e a e ) X 2 (Q 2 )H li2 (x, Q 2 ), 
J- 3 NC (x,Q 2 ) = -2stgn(Q e ){\Q e \(a e + Xv e ) X (Q 2 )G 3 (x,Q 2 ) 



+ 



2v e a e + X (v 2 + a 2 c )] X 2 (Q 2 )H 3 (x, Q 2 )} 



Here structure functions Fj, Gi and Hi describe the hadronic tensor respectively for the 7, 
7 — Z and Z exchange, A = £Qe/\Qe\, £ is the lepton beam polarisation, v e and a e are the 
vector and axial-vector couplings of the lepton to the Z boson: 

v e — 1 — 4|<2e| sin 2 Ow, tie = 1) 

6w is the weak mixing angle, Q e is the lepton charge, Q e = —1, and 

2 G, M 2 Q 2 
X XW ' V2 8na Q 2 + M 2 ' 



z 



with the Fermi constant, G M = 1.16639 ■ l(T 5 GeV~ 2 . 
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The 'radiator' functions <S,- and the function T IR are: 



Si(y,Q 2 ,yh,Ql) 



S3(y,Q 2 ,yh,Ql) 



F IR (y,Q 2 ,y h ,Ql) 
where 

A 2 -- 
B 2 - 

C 2 - 

with 



1 



Ql-Q 2 , (Q 2 + 2m*)(Ql-2mi 



~m 2 (Qf l ~2m 2 



+ 



B 2 



Ql-Q 2 



s <-> -x + 



S2(y,Q 2 ,yh,Ql) = 



[M 2 (Q 2 h + Q 2 )-XS h } 



Q 2 h [S(S - S h ) + X(X + S h ) - 2M\Q{ + 2m 2 )} 



(Ql-Q*)VC~2 
+2m 2 [(S - S h )(X + S h ) + SX] 

-2m 2 -^j- 2 [S(S - S h ) - M 2 Q 2 h ]) - (s <-> -x) - 



2M 2 

7^2 



2QUQl + 2m 2 )(S + X) 



Ql-Q 2 



-2XQ 2 h -S h (Q 2 h + Q 2 ) 



-2m 2 Q 2 h ^(2S - S, 



Q 2 + 2m 



2/ X 



s <-> -xj, 

2 ( B 1 B 2 



Q 2 ~ Ql \VCi \fc~ 2 



m 



3/2 1 £,3/2 



A/ = A u 

2M 2 Q 2 (Q 2 - Ql) + XiStQl - S h Q 2 ) 



+ SQ 2 (S l - S h 



-B^S^-X), 



[XQl - Q 2 (S - S h )] 2 + Am 2 [(S t - S^Ql - S h Q 2 ) 
-M 2 (Q 2 h -Q 2 ) 2 } ee d[S <->-*], 



X = S(l-y) 



2 ME' 

Syh 



and the hadron defined variables, Xh,yh are given by the following equations 

M 2 h = M 2 + Sy h (l-x h ) 
Ql = Sx h y h 
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10 Appendix C 



Below we list the relations between the variables used in the MT (Appendix A) and in the 
D schemes (Appendix B): 



X l = Sf + AM 2 Q 2 (2M|U|) 2 




Mo and Tsai scheme 

E s 

E p 

2Mu 

-g 2 



-x 2 



M 3 



-^^E-E' + Ml Z h o Z ) ) cosfl* 



d (cos9 k ) uj 



q 4 (uq— I u | cos^fe) 
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Figure Captions 

1. Feynman diagrams for the deep inelastic scattering in the one photon exchange ap- 
proximation (a) and the lowest order radiative processes: real photon bremsstrahlung 
from the charged lepton (b), vertex correction (c), vacuum polarisation (d), double 
photon exchange (e), hadron current corrections (f,g). In the MT scheme evaluated 
were diagrams b-d. 

2. Range of kinematical variables from which the radiative tails contribute to the cross 
section measured at the point A(Q 2 , v). 

3. Double photon exchange (a) and hadron current corrections (b,c) in the Dubna scheme. 

4. Infrared cut-off parameter A dependence of the FERRAD35 results obtained for the 
280 GeV muon - proton scattering. Radiative correction factor rj is defined in eq.(f§). 

5. Radiative correction factor r\ calculated in FERRAD35 (open symbols) and TERAD86 
(closed symbols) for the muon - proton scattering at 280 GeV. 
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6. Ratio of the radiative correction factors rj calculated in FERRAD35 (t)f) and TERAD86 
(tjt) for the muon - proton scattering at 280 GeV (closed symbols). The open symbols 
give the r] F /r] T ratio when only e + e~ and contribute to the vacuum polarisation 
(fig. Id) in FERRAD35. 

7. Definition of kinematic variables describing the hard photon emission in MT, cf. Ap- 
pendix A (from JT7]]). 

8. The coordinate system used in the integration over the solid angle of the photon in the 
formula A. 24 in ||, quoted in Appendix A (from ||). 
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